Preprint DFPD/93/TH/45 
June 1993 



Helicity in Classical Electrodynamics and Turbulence 



Aleksej V. Chechkin,Q 

Scientific and Technological Centre of Electrophysics 
Ukrainian Academy of Sciences, P/O Box 9^10, 3rd Valter St., Kharkov, 

310108 the Ukraine 

Dmitrij P. Sorokinfj Q 

Dipartimento di Fisica "G. Galilei", Universita di Padova 

Instituto Nazionale di Fisica Nucleare, Sezione di Padova 
Via Marzolo 8, 35131 Padova, Italia 

and 

Vladimir V. Yanovsky* 

Scientific and Technological Centre of Electrophysics 
Ukrainian Academy of Sciences, P/O Box 9410, 3rd Valter St., Kharkov, 

310108 the Ukraine 

Abstract 



For the electromagnetic fields, hydrodynamic media and turbulent flows we consider the 
relationship between a topological characteristic of vector fields known as helicity and the 
angular momentum of the medium, and discuss, in this respect, the problem of helicity 
and angular momentum transfer from the electromagnetic field to a medium. 
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1 Introduction 



Taking into account topology is one of the most important aspects in studying various 
physical objects. Nontrivial topological properties guaranty the stability of classical and 
quantum states of a physical system, cause the possibility of existing in nature such sta- 
ble point-like objects as vortices, solitons, monopoles and instantons, and may be a cause 
and explanation of several physically observed phenomena (Abricosov vortices, Aharonov- 
Bohm effect, quantum Hall effect etc.) ([1], |2| and references therein). The topological 
characteristics of the system are topological invariants conserved under continuous defor- 
mations of system parameters. From physical point of view these topological invariants 
should correspond to physical observables, i.e. integrals of motion or conserved charges, 
characterizing the dynamics of the system. In this respect establishing the relationship 
between topological and physical characteristics has proved to be important for deeper 
understanding of physical properties and the evolution of the system. 

Among the topological invariants having been engaged in various fields of modern 
physics the Hopf invariant takes one of the first places. In 3-dimensional space it 
characterizes, in particular, the linking and the knot number of the integral lines of a 
vector field A a (x) (a, (3 = 1,2,3) describing a physical system. The helicity of A a (x), 
determined by the integral 

h = — [ d 3 xe a ^A a dBAy (1) 

(where e^ 1 is the unit skew-symmetric tensor with e 123 = 1), has the direct relation to 
the Hopf invariant ||. The integral (1) does not depend on the (local) metric properties 
of the 3-dimensional space and, hence, characterizes global properties of the space and 
the vector field. 

The properties of helical fields have been studied in detail, for example, in the theory 
of magnetic dynamo |4j] and hydrodynamics ||. For instance, it is well known that for 
being able to generate a large-scale magnetic field by a small-scale motion of conducting 
fluid it is necessary for the turbulence to possess helicity. One more example of helicity 
manifestation is the generation of large-scale structures in hydrodynamic turbulence ||. 
If a small-scale turbulence of a fluid or gas possess helicity, then large-scale instability 
leading to the generation of the vortex structures arises. The generation of the large-scale 
structures is caused by the mirror noninvariance of the turbulence. 

In plasmas the problem of current maintenance stimulated great interest in the pos- 
sibility of generating the current by injecting helical electromagnetic waves 0, the prob- 
lem being mainly considered in a way similar to the dynamo effect and examined with 
magneto-hydrodynamic (MHD) equations. 
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In quantum field theory formulated in D = 2 + 1 space-time Eq. (1), called the Chern- 
Simons term (or action), is the basis for developing so called topological quantum field 
theories ||. Being invariant under (abelian) gauge transformations of the vector field 
(A a (x) — > A a (x) + d a (p(x), where <£>(x) is an arbitrary function which at the infinity 
tends quickly enough to zero) the Chern-Simons term allows one in D = 2 + 1 to consider 
massive vector fields without breaking the gauge invariance of the theory || . An approach 
to the description of particles with fractional statistics and spin fLOf , which in their turn 
may occur to be relevant to the problem of high-T c superconductivity JXTJ] , is also based 
on the interaction of matter fields with gauge fields described by the Chern-Simons action 
(1) (see [12 1 and references therein). 

In view of the variety of ways the helicity manifests itself in physics and for deeper 
understanding of this feature it seems of interest to study the meaning of the invariant 
(1) for such well-studied physical objects as electromagnetic fields in 4-dimensional space- 
time, hydrodynamic media and turbulence, and to find the connection of h with observable 
physical characteristics of the system, thus establishing the role of conservation laws 
in helicity transfer from the field to the medium. Since the helicity characterizes the 
'twisting' of the vector fields it is natural to assume that the helicity is connected with 
the angular momentum of the system. For D = 2 + 1 Chern-Simons field theory this is 
well understood, but for ordinary electromagnetic fields and turbulent flows, as we are 
aware, these problems have not been discussed in detail yet. 

In the present paper we consider the relationship between the helicity (1) and the 
angular momentum of the electromagnetic field, hydrodynamic and turbulent flows, and 
discuss the problems of angular momentum and helicity transfer from the electromagnetic 
field to the medium. It is shown (Sections 2, 3) that the trace of the tensor of the spin 
momentum flow density of the electromagnetic field is equivalent to the helicity density of 
the electromagnetic field. Physically it means that helicity density characterizes the value 
of the spin component of the angular momentum flow passing per unit of time through the 
unit surface normal to the direction of electromagnetic field propagation (characterized 
by the Poyting vector) in the observation point. As is well-known, the projection of the 
angular momentum onto the momentum of the electromagnetic field is an invariant of 
the Puancare group called the helicity (not to be confused with (1)!) of the photon. To 
distinguish these two helicities we denote the first of them by h and the second one by s. 
In classical electrodynamics s-helicity corresponds to the degree of electromagnetic wave 
polarization. As we shall see, Eq. (1) determines the degree of the circular polarization 
of the electromagnetic wave and is proportional to the mean value of angular momentum 
projection on to the Poynting vector. In the case of the plane wave Eq.(l) indeed coincides 



2 



with the expression for the classical analogue of the photon helicity s, and, therefore, is 
conserved (physically observable) quantity 

In Section 4 by analogy with the electromagnetic case we consider the relationship 
between the helicity and angular momentum for turbulent flows, and introduce the kinetic 
definition of the helicity of particle systems. This kinetic definition is similar to the 
hydrodynamic one, but it allows us to get deeper insight into the microscopic nature of 
the /i-helicity. 

In Section 5 a problem of angular momentum and helicity transfer from the field to 
the medium is discussed by use of a simple example of a charged particle propagating in 
a circular polarized plane electromagnetic wave. It is shown that the particle can acquire 
angular momentum from electromagnetic waves with nonzero helicity. 



2 The angular momentum of the electromagnetic field 

A density tensor of the angular momentum of the electromagnetic field may be constructed 
by the analogy with the orbital momentum tensor of the classical particle ]T3| : 

= Xi T* - Xj T*, (2) 

where 

Ta = ^ c ( F ?F kj - \ 9lJ F kl F kl ) (3) 

is the symmetric energy-momentum tensor and Fik = d{A^ — diAk is the electromagnetic 
field stress tensor, is the Minkowski metric with a signature chosen to be (+, — , — , — ), 
and i,j, k, ... = 0, 1, 2, 3. 

However, more natural and favorable (for our consideration) is to write down a density 
tensor of the angular momentum from fundamental principles of the symmetry of elec- 
trodynamics under the Lorentz group transformations i.e. by varying the electromagnetic 
field Lagrangian in compliance with the Noether theorem (see, for example, [0). In this 
case the angular momentum density tensor has the form 

*4 N )k = [M* + ^(x^A^F* - xj&AAF 1 *)] + ^(-^ fc + AjFf), (4) 

where the first two terms correspond to the density of the orbital momentum and the 
last two terms describe the spin momentum of the electromagnetic field. It is just the 
manifest representation of the angular momentum as the sum of the orbital and the spin 
part enables one, as we will see below, to understand the physical meaning of Eq.(l). 
These two definitions, (||]) and (§]), of the angular momentum differ by a total derivative 

M\f )k - M% = -L^x^F 1 * - x,AF lk ), (5) 
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and since physical observables are spatial integrals of the conserved current densities and 
under the condition that the fields tend quickly enough to zero at infinity, the integrals 
of (HD and (JD taken over the volume of 3-dimensional space coincide. However, if the 
fields tend to zero rather slow or spread to infinity (as, for example, the plane waves), 
then taking into account boundary conditions becomes essential for the computation of 
the integrals of @ and (£|). In overwhelming majority of cases this leads to interesting 
physical observations concerning (global) properties of the physical system. On the other 
hand the neglect of boundary conditions may result in incorrect physical conclusions. 

The simplest example is the calculation of the projection of the angular momentum on 
to the momentum of a plane electromagnetic wave. If for this calculation we use Eq. (^), 
then the result will be identically zero. Indeed, the Poynting vector of the plane wave has 
the form 

P a = [ d'xT° a (x) = d*x(E 2 + H 2 ), (6) 

J Q71C J 

where n a is the unit vector along the direction of wave propagation, and E a = F a0 and 
H a = ^e al 3-yF al3 are electric and magnetic field strength, respectively, satisfying (for the 
plane wave) the relation 

H = n x E. (7) 
The angular momentum calculated with the Eq. (0) has the form 

M a = -—e aPl n p f d 3 xx 7 (E 2 + H 2 ), (8) 

87TC J 

Thus, the helicity 

MP 

s = w {) 

of any plane electromagnetic wave identically terns to zero if one substitutes into @ 
Eq. (P) and (||), which runs counter the observation: in general the plane wave has 
nonzero circular polarization whose degree is characterized just by the magnitude of s 
(or, another words, it is well-known that the photon has nonzero helicity). 

On the other hand the computation of s with the use of Eq. (0) of the angular momen- 
tum density tensor gives the correct result and allows one to find the connection between 
h and s helicity of the electromagnetic wave. 

It is convenient to carry out the computation in the Coulomb gauge 

A = 0, ^A«(x) = 0, (10) 

then the only contribution into the value of s (§) is given by the spinning part of the 
angular momentum (|j), and Eq. @ (with taking into account (^)) takes the form 

s = —e al3y [ d 3 xn a A g E~ = — f d 3 xA a H a = — / d 3 xA(V x A), (11) 

4ttc 11 J Attc J 4ircJ ' v ' 
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which coincides with the /i-helicity (Eq. (1)) of the electromagnetic field. 

Thus, the s-helicity (or circular polarization) of the plane electromagnetic wave is 
determined by Eq. (1), which allows one to make a conclusion, that the integral lines 
of electromagnetic field with nonzero circular polarization have the nontrivial topological 
structure characterized by a linking number. 



3 Angular momentum flow and h— helicity 

Let us show that the /i-helicity of an arbitrary electromagnetic field is equal to the trace 
of the spin component of the angular momentum flow. 

It follows from the conservation law for the angular momentum (||) 

WS?o - d J M iZ = ( 12 ) 
that the density of the spin momentum flow has the form 

G a p = —~ Za-fSA^Fsp, (13) 

Z7TC 

and 

G = TrG aP = -—e ajS A 7 F 50 = -^A(V x A). (14) 

Z7TC 47TC 

Thus, the density of the spin momentum flow along the direction of electromagnetic field 
propagation (i.e. the density of s-helicity flow) is determined by the /i-helicity density, 
and Eq. (1) characterizes the average value, over spatial volume, of the s-helicity flow 
density. The topologically nontrivial configurations of the electromagnetic field manifest 
themselves through the nonzero helicity values. 



4 Helicity in turbulence 

Let us now proceed with studying the relationship between the helicity and angular mo- 
mentum of turbulent flows. As we have already mentioned in the Introduction the notion 
of the helicity of random fields plays an important role in the theory of turbulence; thus, 
establishing the connection between these two physical notions may occur to be useful for 
better understanding the turbulence phenomenon. 

We shall consider this problem by the example of the most simple type of turbulent 
motion being homogeneous in space and time and isotropic in space. Just that every case 
has been considered for studying problems of magnetic field generation (a-effect) and 
the generation of large scale structures by small scale helical fields . 
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The pair correlation tensor of the velocity field V(r, t) 

(V i {r,t)V j (r f ,f))=B ij (r,T>,t,lf). (15) 

('( ) 'denotes statistical averaging) is the most important characteristic of turbulence. So 
let us discuss this characteristics in more detail. The requirement for the tensor (|15|) 
to conserve the form under rotations and shifts in the space leads to the following most 
general form of the correlation tensor [ITH 



where R = r — r' and r = t — t', A(R, r) and B(R, r) are scalar functions and G(R, r) 
is a pseudoscalar function. All this functions are expressed through the corresponding 
correlation functions of V(r, t). For example, 

4(0,0) = ^(V i (T,t)V i (T,t)) 

has the physical meaning of the average energy density of the fluid, and 

G(0,0)4(V(r,i)(VxV(r 1 t))) (17) 
o 

determines the helicity density of V. 

The physical meaning of B(R, r) is more complicated and we will not discuss it here. 

One can easily see that the form of the right hand side of (|1^) is conserved not only 
under the spatial rotations and shifts but under the reflections as well. But since the 
statement that turbulence is invariant under some group transformations means that all 
correlation characteristics of turbulence (and not only ,for example, (0)) are forminvari- 
ant under this transformations, and ( |T7j ) is indeed reflection noninvariant, it follows that 
turbulence with nonzero G(0, 0) is reflection noninvariant. 

Let us show now that in agreement with the results of the previous sections concerning 
the helicity of the electromagnetic field the helicity of the turbulent velocity field also has 
(up to a dimensional factor) the "dynamical" meaning of the projection of the total orbital 
momentum onto the total momentum of the turbulent medium. 

Assuming for simplicity that the density p of the fluid is constant we wright down the 
projection in the form 

* = y</ dr[r x V(r,t)] J dr'V(r', t)), (18) 

where V is the volume of the medium. Then, using (|16D , the equation fll8|) may be reduced 
to the form 

s = p J dRR 2 G(R, 0) = p J dRR 2 (V(R, 0)V x V(R, 0)) . (19) 
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If for R larger then a correlation radius G(R, 0) tends to zero quick enough we get from 



(|D that 

s~(V(r,t)VxV(r,t)), (20) 

Eqs. (0), ( |2"0"D represent the relationship we have been looking for. The observation 
of this fact points to the close connection between topological properties and dynamical 
characteristics of turbulence. Such a connection plays an important role for understand- 
ing the processes of helicity (or angular momentum) exchange between turbulent flows 
and external fields acting on the former, the generation of various turbulence structures 
possessing helicity, etc. [|], 0. A simple example of helicity and angular momentum 
transfer is considered in the next section. 

Using the results obtained we may go even further and introduce the notion of helicity 
for a system of particles by defining a microscopic helicity as 



s 



M 



£[r a x P Q ]£P6, (21) 



a b 

where r a and p a are the radius vectors and the momentums of the particles, a, b = 1, N 
and iV is the number of particles. Using the microscopic phase density |16[ 



Af (r, p, t) = £ S (r - r(t)) 6 (p - p(t)) (22) 

a 

we can rewrite fl2"T|) as 

~ J cMcM'ip x p'](r - r')Ar(r,p,t)AT(r',p',t), (23) 



s M 



where dQ is the measure of integration over the phase space. 

Taking into account the fluctuation of the microscopic phase density 

5J\f = A/" - AT, 

where A" = nf(r,p,t), f(r, p, t) is the one particle distribution function and n = M /V 
is the average particle density in the system, we write down the average value of the 
microscopic helicity as follows: 

s = (s M ) = \J cMcM'ip x p'](r - r')/(r, p, t)f(r>, p', t) 

+ ^ J dttdn'lp x p'](r - r') (8X(T,p,t)8X(r f ,p',t)) . (24) 

The first term in (^) describes a "hydrodynamic" part of the average helicity and can 
be rewritten in the form (compare with (|18|) ) 

n 2 m 2 J drdr'[V(r,t) x V(r, t)](r - r'), (25) 
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where V(r, t) — — J dppf(r, p, t) and m is the particle mass. 

The second term characterizes a "correlation" part of the average helicity. If the 
correlation function of the phase density fluctuation is isotropic the correlation part of the 
average helicity turns to zero. One may also show that the correlation function of a weak 
interacting homogeneous classical gas does not contribute, in the Landau approximation, 
into the helicity, and, hence, an inhomogeneous state should be taken into account. Thus, 
for these cases the kinematic definition of the helicity of a particle system introduced 
above coincides with the hydrodynamic definition of the helicity. In general case the 
correlation part of the helicity may cause new correlation effects, but this point requires 
additional studies. 



5 Angular momentum transfer from field to medium 

In previous sections we have elucidated the connection between the angular momentum 
and the /i-hehcity of the electromagnetic field and medium. 

In many physical applications a problem arises to determine the value of energy, mo- 
mentum and angular momentum transferred from the field to the medium. This problem 
is important for understanding and predicting the behavior of the system. 

However, while, for example, in plasmas and turbulence the energy and momentum 
transfer have been studied in detail |[L7|| , the problem of angular momentum transfer has 
not been developed in the same extent, though it has been under consideration since 
Sadovski found the effect of turning a small plate by an electromagnetic wave [18| . 



We see that the problem of angular momentum transfer is related to the problem of 
the helicity transfer from fields to the medium discussed, for example, in [^]. 

The problem of angular momentum transfer from electromagnetic fields to particles 
and consequences of this effect for the electromagnetic pumping up of the angular mo- 
mentum into plasmas were discussed, in particular, in Ref. ||19|| . Here we consider this 
question from somewhat different point of view based on the relationship between the 
angular momentum and the /i-helicity, and, in contrast to [Hj, demonstrate that the 
charged particles may acquire orbital momentum from plane electromagnetic waves with 
nonzero helicity and that this effect may result in pumping up the angular momentum 
into plasma. 

We restrict ourselves by the consideration of the simplest possible model of a system of 
non-interacting charged particles, so that the problem is actually reduced to considering 
the motion of a single particle. 

Let a particle with a charge e and a mass m be placed into the field of a circular 
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polarized plane wave propagating along the axis x and having electric strength 

E = e y E cosu;(£ ) + ae z E sinu;(t ), (26) 

where Eq is the wave amplitude, u is the electromagnetic field frequency and a = ±1 
denotes left or right polarization. 

Note that due to Eqs. @, (|ll|) the angular momentum of the plane wave is directed 
along the wave momentum (axis x), and its value coincides with the helicity of the elec- 
tromagnetic wave, the helicity density of the latter being 

„ a . _ , aEl ae 

G = - — AV x A = — ± = — , (27) 
Aire Aituj uj 

E 2 

where e = is the energy density of the wave. If we considered a wave with arbitrary 
polarization the value of a in (p?D would have determined the corresponding polarization 
properties of the wave. For instance, if the wave is linear polarized, a = 0, and, hence 
h = s = 0, which testifies to the triviality of the corresponding topological characteristics 
of the linear polarized wave. 

The equation of motion of the particle in the electromagnetic wave is nonlinear, and, 
in general, finding an exact solution of the equation is nontrivial problem.[| In view of this, 
and since here we do not aim at carrying out the rigorous derivation of the expression 
for the particle orbital momentum acquired, we shall consider the process of angular 
momentum transfer from the wave to the particle using semiqualitative consideration 
analogous to that used, for example, in the computation of the strength of high frequency 
pressure . 



We look for the solution as a series expansion in the field amplitude and restrict 
ourselves to the nonrelativistic approximation. Particle motion consists in quick rotation 
in (y, z) plane and slow motion averaged over the oscillations. Let us further assume that 
the particle velocity (v) averaged over the oscillations is directed along wave propagation 
(axis x). We shall see below that for particle to get an angular momentum from the wave 
it is crucial that the former accelerates along the wave. Slow variation of (v) with time 
may be caused, in the first place, by the momentum transfer from the scattered wave to 
the particle, or by a weak constant electric field. 

Then from the equation 

m^ = eE(<r),f) + ^v)xH(<r),t) (28) 



-'^Note that the well known particular solution of the equation, which describes particle motion along 
a circle perpendicular to the wave momentum, implies that the particle had the angular momentum a 
priori, and the latter does not vary in the course of motion. Hence, this solution is not one we are looking 
for. 
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for the oscillating part f of the particle radius vector we get 

f = o 6E ° (ri^ Mo cos (1 - M)arf + a e z E sin (1 - ^-)ut). (29) 

mu! 2 (l — — ) c c 

««> = |<v>|) 

Eq. ( |29| ) may be directly obtained by solving Eq. ( p8|) in the coordinate system where 
(v) is equal to zero in a given time moment and by passing to the laboratory system 
with making use of the Lorentz transformation (Doppler effect) in the nonrelativistic 
approximation. 

Momentum variation with time is expressed by the formula 

dM . d 2 r, , . 

if = |r x m ^ (30) 

Taking into account slow dependence of (v) on time we get from (^) and Q30D that 

dM = An^d{vl G 
dt mew 2 dt 

We see that the rate of momentum variation is proportional to the helicity density 
G and particle's slow acceleration along the wave. This testifies to the existence of the 
effect of momentum transfer from an electromagnetic wave to a charged particle, and 
indicates that such a transfer is possible if only the wave possesses the helicity, the process 
being governed by the angular momentum conservation law. The simple example having 
been considered shows a way the topologically nontrivial fields affect the dynamics of the 
charged particles. 

For plasmas this effect means the possibility of pumping up the angular momentum 
into the plasmas by the helical fields. The angular momentum of the plasma will increase 
until the particle acceleration along the wave is compensated by the friction force, and 
until the increase is terminated because of the angular momentum loss caused by particle 
radiation. 



6 Conclusion 

For electromagnetic fields, turbulent media and particle systems we have established the 
relationship between the /i-helicity (1), being a characteristics of topologically nontrivial 
vector fields, and the s-helicity (|^) which characterizes the value of the projection of the 
angular momentum onto the momentum of the medium. This testifies to the close connec- 
tion between corresponding topological and dynamical properties of the physical objects, 
which should be taken into account when studying the effects of angular momentum and 
helicity exchange between fields and media. 
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For hydrodynamics and turbulence the momentum transfer between fields and media 
points to the "intertwining" of the helicity of the field with that of the medium velocity 
field, which may point to the existence of a conservation law for a topological charge 
characterizing the field-medium system CIS db whole. 

Note added. When this work was completed the authors became aware of Ref. [ET[|, 
where the connection between the two types of helicity of the electromagnetic field was 
pointed out in the framework of a topological electromagnetism. 
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